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ABSTRACT 

Various resting-state fMRI (R-fMRI) measures have been developed to characterize intrinsic 
brain activity. While each of these measures has gained a growing presence in the literature, 
questions remain regarding the common and unique aspects these indices capture. The present 
work provided a comprehensive examination of inter-individual variation and intra-individual 
temporal variation for commonly used measures, including fractional amplitude of low frequency 
fluctuations, regional homogeneity, voxel-mirrored homotopic connectivity, network centrality 
and global signal correlation. Regardless of whether examining intra-individual or inter-individual 
variation, we found that these definitionally distinct R-fMRI indices tend to exhibit a relatively 
high degree of covariation. When taken as a measure of intrinsic brain function, inter-individual 
differences in concordance for R-fMRI indices appeared to be stable, and negatively related to 
age (i.e., functional concordance among indices decreases with age). To understand the 
functional significance of concordance, we noted that higher concordance was generally 
associated with higher strengths of R-fMRI indices, regardless of whether looking through the 
lens of inter-individual (i.e., high vs. low concordance participants) or intra-individual (i.e., high 
vs. low concordance states identified via temporal dynamic analyses) differences. Finally, 
temporal dynamics analyses also revealed that high concordance states are characterized by 
increased within- and between-network functional connectivity, suggesting more general 
variations in network integration and segregation. The current study draws attention to 
questions regarding how to select an R-fMRI index for usage in a given study, as well as how to 
compare findings across studies that examine inter-individual or group differences using 
different indices. Additionally, our work suggests global neural signals exist in the brain, and 
their spontaneous variations over time result in fluctuations in the connectedness of brain 
regions. 
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INTRODUCTION 

Intrinsic brain activity (IBA) is posited to play a central role in human brain function. In particular, 
it is believed to help maintain a state of preparedness through the continuous representation of 
high probability patterns of functional interactions at the local and systems levels, regardless of 
stimulation level or task demands (Raichle, 2010, 2015). In support of the proposed functional 
significance of IBA, prior work has linked intra-individual variations in behavior to the state of 
spontaneous activity in task-relevant systems, as well as those that can interfere (Fox et al., 
2007; Kelly et al., 2008). Additionally, a growing literature has worked to link inter-individual 
variations in trait IBA properties (e.g., functional connectivity) to differences in complex 
behavioral phenomena (e.g., personality, social reciprocity, intelligence, anxiety) among 
individuals (see reviews in Di Martino et al., 2014a; Fox and Raichle, 2007; Kelly et al., 2012). 
As efforts to characterize IBA differences among individuals and clinical groups has increased, 
a growing number of measures have been proposed - each posited to reveal a distinct aspect 
of intrinsic brain function (see reviews in Craddock et al., 2013; Margulies et al., 2010; Zuo and 
Xing, 2014). While definitional differences exist among these various measures from a 
theoretical perspective, practical distinctions and commonalities among these measures are yet 


to be fully explored empirically. 


The most fundamental IBA measures are those focused on regional fluctuation properties. 
Among them, amplitude-based measures have gained particularly attention and can be carried 
out in either the temporal [e.g., resting state fluctuation amplitude (RSFA, Kannurpatti and 
Biswal, 2008)] or frequency domains [e.g., amplitude of low frequency fluctuations (ALFF, Zang 
et al., 2007); fractional ALFF (fALFF, Zou et al., 2008); Hearst Component (Maxim et al., 2005)]. 
Complementing these regional measures, is functional connectivity, which indexes the inter- 
regional synchronization of IBA using measures of temporal correlation or spatiotemporal 


component decompositions, such as independent component analysis (ICA, Beckmann et al., 


2005)]. Building on top of the inter-regional synchronization measures, are those focused on 
particular theoretical constructs, such as local synchronization [e.g., regional homogeneity 
(ReHo, Zang et al., 2004; Zuo et al., 2013)], interhemispheric synchronization [e.g., voxel- 
mirrored homotopic connectivity (Anderson et al., 2011; VMHC, Zuo et al., 2010b)] or network 
centrality [e.g., degree centrality (DC, Buckner et al., 2009; Zuo et al., 2012)]. While each of 
these measures is widely employed to examine individual and group differences, they are 
commonly studied in isolation of one another. As a result, questions remain regarding how 
similar or distinct the various IBA indices are, especially in the absence of an integrative 


neurophysiological perspective of their underlying mechanisms. 


Realizing the potential for differences and overlaps among these indices, some studies have 
attempted to include a range of measures in their examination of individual- or group- 
differences. For example, Yan et al. (2013a; 2013b) examined inter-individual differences in 
ALFF, fALFF, ReHo, DC, intrinsic functional connectivity (iFC) with a focus on devising 
strategies for reducing and standardizing their variability. Chen and Xu et al (2015) 
systematically investigated various IBA metrics regarding their intra-individual, inter-individual 
variability and test-retest reliability using richly sampled R-fMRI datasets. In a recent study, 
fALFF, ReHo and DC were observed to converge at the left angular gyrus regarding their 
across-individual correlation with the working memory maintenance (Yang et al., 2015). 
Similarly, an array of IBA measures was reported to converge their alterations at insula and 
posterior cingulate cortex in patients with Autism Spectrum Disorders (Di Martino et al., 2014b). 
Such convergence in findings may suggest the presence of a more fundamental mechanism 


underlying the different IBA measures and their variations across individuals. 


The goal of the present work is to provide a comprehensive understanding of interdependencies 
among different IBA measures within and across individuals. In service of this goal, we identified 
three potential strategies: 1) Examine how concordant differing indices are with respect to their 
variation across voxels. From this perspective, prior work has found that ALFF is strongly 
coupled with fALFF (Zou et al., 2008; Zuo et al., 2010a), and also positively correlated with 
ReHo across voxels (Nugent et al., 2015; Yuan et al., 2013). In a recent study, Aiello et al. 
(2015) examined the inter-correlation between fALFF, ReHo and DC across voxels, finding high 
pair-wise correlations between the 3 measures across voxels. 2) Examine how concordant 
different indices are with respect to their variation from one individual to the next (i.e., are they 
capturing unique inter-individual variation vs. shared). Aiello et al. (2015) also provided initial 
insights in this regard, as they comprehensively examined correlations between fALFF, ReHo 
and DC across participants, finding high correlation for all the gray matter voxels. Beyond that, 
Yang et al. (2016) used a fully data-driven discovery-replication framework to examine inter- 
individual differences in 6 IBA measures (ALFF, fALFF, ReHo, ReHo-2, DC and eigenvector 
centrality) and 6 morphology measures, mapping a set of interdependencies of the IBA metrics 
and their morphological associations. 3) Examine how concordant differing indices are with 
respect to their variation over time (e.g., from one time window to the next within a given 


individual), this perspective is yet to be investigated. 


While optimal methodologies are still being determined (Hindriks et al., 2016), the temporal 
dynamic perspective of IBA measures makes it possible to determine how differing measures 
can be coupled together within an individual. A growing number of studies have demonstrated 
temporally dynamic changes in iFC patterns (Allen et al., 2014; Chang and Glover, 2010; 
Handwerker et al., 2012; Hutchison et al., 2013; Kang et al., 2011; Kiviniemi et al., 2011; Smith 
et al., 2012; Yang et al., 2014)). Based on the spatiotemporal structures of iFC, brain activity 


can be decoded into distinct brain states (Allen et al., 2014; Yang et al., 2014; Zalesky et al., 


2014). Dynamic functional network analysis also showed that the dynamic FC fluctuations 
coincide with periods of high and low network modularity in IBA (Betzel et al., 2016). However, 


the temporal dynamics of regional IBA as well as their interdependencies remain unexplored. 


The present work made use of the publicly available Enhanced Nathan Kline Institute - 
Rockland Sample dataset, which provides resting state fMRI data with a relatively high temporal 
resolution (TR = 0.645s) compared to traditional imaging sequences (thereby making it more 
optimal for temporal dynamics anaylses). Similar to prior work, we first examined the inter- 
correlation and concordance of differing R-fMRI measures (ALFF/fALFF, ReHo,DC and УМНС) 
across participants. We then investigated the temporal dynamics of these measures and the 
extent to which their variations over time are coupled with one another. To determine if a 
common driving force (either physiological or neuronal) drives their temporal dynamics, we also 
investigated the time varying pattern of global signal correlation (GSCorr). We also tested for 
any potential relationships with physiological noises (respirational volume and heart rate). 
Finally, to examine their potential value in neuronal applications, we investigated how coupling 


relationships change with age across the lifespan. 


MATERIALS AND METHODS 

Participants and Imaging Protocols 

Multiband fast-sampling R-fMRI data (TR = 0.645s, 10 minutes/900 volumes) of 173 
neurotypical individuals ages between ages 8.0 and 86.0 with quality pass datasets (mean age: 
44.5; 117 females) were selected from the publicly available Enhanced Nathan Kline Institute - 
Rockland Sample (Available via International Neuroimaging Data-sharing Initiative: 
http://fcon_1000.projects.nitrc.org). Images were acquired in а ЗТ Siemens TIM Trio scanner 


using the following multiband echo-planar imaging (EPI) sequence: TE = 30 ms, flip angle = 60*, 


slice thickness = 3.0 mm, field of view = 222mm, matrix size = 74 x 74,TR = 645ms; no gap, 
resolution = 3.0 x 3.0 x 3.0 тт?. During the data acquisition, participants were instructed to 
simply rest with eyes open (Nooner et al., 2012). For spatial normalization and localization, a 
high-resolution T1-weighted magnetization prepared gradient echo image (MPRAGE) was also 
obtained for each participant. The Nathan Kline Institute institutional review board approved the 


submission of anonymized data obtained with written informed consent from each participant. 


Preprocessing 

Unless otherwise stated, all preprocessing was performed using the Data Processing Assistant 
for Resting-State fMRI (DPARSF, Yan and Zang, 2010, http://rfmri.org/DPARSF), which is 
based on Statistical Parametric Mapping (SPM8) (http://www.fil.ion.ucl.ac.uK/spm) and the 
toolbox for Data Processing & Analysis of Brain Imaging (DPABI, Yan et al., in press, 
http://rfmri.org/DPABI). The first 10 s of data were removed to allow data to reach equilibrium, 
leaving a total of 884 volumes for final analysis. The functional volumes for each subject were 
motion-corrected using a six-parameter (rigid body) linear transformation with a two-pass 
procedure (registered to the first image and then registered to the mean of the images after the 
first realignment). Individual structural images (T1-weighted MPRAGE) were co-registered to the 
mean functional image after realignment using a 6 degrees-of-freedom linear transformation 
without re-sampling. The transformed structural images were then segmented into gray matter 
(GM), white matter (WM) and cerebrospinal fluid (CSF) (Ashburner and Friston, 2005). The 
Diffeomorphic Anatomical Registration Through Exponentiated Lie algebra (DARTEL) tool 
(Ashburner, 2007) was used to compute transformations from individual native space to MNI 


space. 


Nuisance Regression 


Recent work has demonstrated that micro head movements, as small as 0.1mm between time 


points, can introduce systematic artifactual inter-individual and inter-group variability in R-fMRI 
measures (Power et al., 2012a, b; Satterthwaite et al., 2012; Van Dijk et al., 2012). Here we 
utilized the Friston 24-parameter model (Friston et al., 1996) to regress out head motion effects 
from the realigned data (i.e., 6 head motion parameters, 6 head motion parameters one time 
point before, and the 12 corresponding squared items) based on recent reports that higher- 
order models demonstrate benefits in removing head motion effects (Satterthwaite et al., 2013; 
Yan et al., 2013a). Head motion was also controlled at the group-level by taking mean 


framewise displacement as a covariate (Satterthwaite et al., 2013; Yan et al., 2013a). 


The signals from WM and CSF were regressed out to reduce respiratory and cardiac effects. In 
addition, linear and quadratic trends were also included as regressors since the blood oxygen 
level dependent (BOLD) signal demonstrates low-frequency drifts. Global signal regression 
(GSR) is not performed because of its controversy in increasing negative correlations (Murphy 
et al., 2009; Weissenbacher et al., 2009) and may distort group differences in iFC (Gotts et al., 
2013; Saad et al., 2012). Temporal filtering was not performed to keep the temporal dynamics of 
time series. After nuisance regression, the functional data were registered into MNI152 space 
with 3mm? cubic voxels by using transformation information acquired from the previous 


DARTEL step. 


R-fMRI Indices for Regional Characteristics and Functional Synchrony 

We examined the interdependency among the following five R-fMRI based indices of intrinsic 
brain activity: 

1) Amplitude measures: ALFF (Zang et al., 2007) and fALFF (Zou et al., 2008). ALFF is the sum 
of amplitudes within a specific low frequency range (0.01 - 0.1 Hz in the current study) from a 
Fourier decomposition of the time course. fALFF is the ratio of the ALFF of a given low 


frequency band (here, 0.01 - 0.1Hz) to the sum of Fourier amplitudes across the entire 


frequency range. ALFF is proportional to the strength or intensity of low frequency oscillations, 
while fALFF represents the relative contribution of specific oscillations to the whole detectable 
frequency range. Given the high colinearity between fALFF and ALFF, we chose to include only 
one of the two measures in the concordance indices. fALFF was chosen because it’s less 
susceptible to artifactual contributions of motion and pulsatile effects as compared with ALFF 
(Yan et al., 2013a; Zou et al., 2008; Zuo et al., 2010a). 

2) Regional homogeneity (ReHo, Zang et al., 2004). ReHo assesses the degree of regional 
synchronization/coherence among fMRI time courses. It is defined as the Kendall’s coefficient of 
concordance (or Kendall's W, Kendall and Gibbons, 1990) between the time series of a given 
voxel with those of its nearest neighbors (26 in the current study). 

3) Homotopic interhemispheric connecitivty. VMHC (Anderson et al., 2011; VMHC, Zuo et al., 
2010b) corresponds to the functional connectivity between any pair of symmetric inter- 
hemispheric voxels - that is, the Pearson's correlation coefficient between the time series of 
each voxel and that of its symmetrical inter-hemispheric counterpart. The resultant VMHC 
values were Fisher-Z transformed. For better correspondence between symmetric voxels, 
VMHC requires that individual functional data are first registered in MNI space and smoothed 
(4.5 mm FWHM) and then registered to a symmetric template. The group averaged symmetric 
template was created by first computing a mean normalized T1 image across participants, and 
then this image was averaged with its left-right mirrored version (Zuo et al., 2010b). 

4) Network degree centrality (Buckner et al., 2009; Zuo et al., 2012). Degree centrality is the 
number or sum of weights of significant connections for a voxel. Here, we calculated the 
weighted sum of positive correlations by requiring each connection's statistical significance to 
exceed a threshold of r » 0.25 (Buckner et al., 2009). 

5) Global Signal correlation (GSCorr): GSCorr is the Pearson correlation coefficient between the 
global mean time series (averaged across all voxels inside the group mask) and each voxel in 


the group mask. These correlation values were then Fisher-Z transformed. 


Dynamic R-fMRI Indices 

Dynamic R-fMRI indices were generated using sliding time-window analysis. First, we applied 
hamming windows (length of 100 TRs, overlapping of 3 TRs)' to BOLD signals to obtain 
windowed time series; second, within each window, we calculated the above-mentioned R-fMRI 
indices (i.e., fALFF, ReHo, VMHC, DC and GSCorr). To characterize the dynamic R-fMRI 
indices, we compute the mean and standard deviation (SD) map across time windows for each 
index. To quantify temporal variations of these indices, we also computed the coefficient of 
variation (CV: mean/SD) map over time for each index. The mean, SD, and CV maps were then 
Z-standardized relative to its own mean and SD across all the voxels within the group mask. 
Across participants, we performed one-sample t-tests on the Z maps of mean, SD, and CV to 


investigate the overall pattern of temporal variation. 


Correlation between Global Mean of R-fMRI Indices 

For static analysis, the global mean of each index was computed within the group mask for each 
participant. Pearson’s correlation coefficient was then computed between pairs of the global 
mean of R-fMRI indices across participants. For dynamic analysis, we calculated the global 
mean of a given R-fMRI measure for each time window. This resulted in one global mean time 
series per index per participant. Within each participant, we calculated the correlations between 
pairs of the global mean time series of R-fMRI indices. These correlations were then averaged 
across participants to provide an estimation of the overall dynamic correlations. To examine if 
the temporal dynamics of these measures are being driven by physiological noises, we 
examined their relationship with each: motion, global signal, respiration volume per time (RVT) 
and heart rate (HR). RVT and HR were extracted by applying AFNI's RetroTS program on the 


physiological data. 


1 Window size of 64 TRs and 128 TRs were also tested, resulted in similar results. 


Voxel-wise Concordance Index 

We calculated the Kendall's W of the five static R-fMRI indices across participants for each 
voxel as a static voxel-wise concordance index. To compare the within-individual dynamic 
concordance to the static functional concordance across individuals, we calculated Kendall's W 
for the five R-fMRI indices across time windows as the dynamic voxel-wise concordance index. 
A one-sample t-test was performed on the Z-standardized dynamic voxel-wise concordance 
index to examine the pattern among subjects. The rank order of the static voxel-wise 
concordance, as well as the rank order of the group t map of dynamic voxel-wise concordance 


index, was calculated to compare their consistence. 


Volume-wise Concordance Index 

For static analysis, we calculated Kendall’s W of the R-fMRI indices across all brain voxels for 
each participant. To investigate the properties of the functional concordance, we compared 
subjects with the highest 20% static volume-wise concordance to those with lowest 20% static 
volume-wise concordance. For each R-fMRI measure, two-sample t-tests were performed 
between the highest 20% and lowest 20% subjects to capture the group differences. For 
dynamic analysis, dynamic volume-wise concordance index was computed the same way as 
static analysis, but for each time window. The mean dynamic spatial functional coupling index 
(across time) was then compared with the static volume-wise concordance index. To 
characterize the properties of high concordance states in contrast with low concordance states 
at the within-individual level, we compared properties of top 20% windows with highest dynamic 
volume-wise concordance to those windows with the lowest 20% dynamic volume-wise 
concordance. For each individual, a mean map of a given measure was calculated for the top 
20% windows, as well as a mean map for the lowest 20% windows. Then paired t-test was 


performed to compare the property difference between the high windows and the low windows. 


Age Effects on Functional Concordance 

Finally, to examine potential role of functional concordance in development, we linked the 
functional concordance indices to age across the lifespan. A general linear model with sex and 
motion as covariates was utilized with DPABI. Specifically, for voxel-wise or volume-wise 


concordance,index, the following regression model was constructed: 


A given measure = b0 + b1xAge + b2xSex + b3xmeanFD + error 


Mean FD was included to control for the residual effect of head motion. We opted to employ 
group-level corrections over scrubbing approaches for the following reasons: (1) recent work 
suggested that scrubbing offers little advantage over group-level motion correction combined 
with Friston 24 model at the individual level. Scrubbing can also be less conservative when 
motion correlates with a between-subject variable of interest, and can be incomplete in their 
corrective effects (Yan et al., 2013a); (2) scrubbing of non-contiguous time points alters the 
underlying temporal structure of the data, which can unintentionally bias estimates of an 
individual's intrinsic activity patterns towards one state over another in temporal dynamics 
studies. The results were corrected for multiple comparisons using Gaussian random field 
theory. In appreciating recent reports on enhancing the primary threshold (P value at individual 
voxel level) (Eklund et al., 2015; Woo et al., 2014), we have set voxel-level P<0.001 (Z>3.1), 


cluster-level P<0.05 for multiple comparison correction. 


RESULTS 
Static and Dynamic R-fMRI Indices 
The first step in the present work was to characterize the temporal variation in our R-fMRI 


indices of interest. In this regard, for each of the indices, we carried out sliding window analyses 


(window size: 100 TRs = 64.5 sec) and then calculated the voxel-wise coefficient of variation 
across windows for each participant. The resultant coefficient of variation maps are depicted in 
Figure 1, along with the traditional static R-fMRI index maps. As expected, the temporal mean 
and static index map were nearly equivalent. Inspection of the CV maps revealed different 
patterns of variation across R-fMRI indices. fALFF appeared to show lower temporal variation in 
high level cognitive regions (e.g., default mode network components, dorsolateral prefrontal 
cortex), and higher variation in lower-level primary sensory and motor areas, as well as 
subcortical regions. In contrast, for ReHo and DC, the variation is high among all the cortices, 
including the default mode network. Of note, CV is not applicable for VMHC and GSCorr as they 


contain negative correlation values. 


Evaluating Concordance among R-fMRI Indices: Global-Level Analyses 

From the perspective of inter-individual variation, all measures were closely coupled with each 
other at the global level. Static analyses that calculated for each participant the mean grey 
matter value for each of the measures of interest, found high correlation between each pair of 
R-fMRI indices across individuals (Figure 2A). Dynamic analyses capitalized on the ability to 
perform within-individual analyses; specifically, at each time window, the mean gray matter 
value for each of the measures was calculated for each participant; a high degree of correlation 
was noted between all pairs of R-fMRI indices over time (mean correlation scores across 
participants depicted in figure) (Figure 2B), indicating common fluctuation patterns underlying all 
the functional aspects. Of note, GSCorr is less correlated with fALFF (mean r = 0.54), but highly 
correlated with DC (mean r = 0.87). Nuisance parameters, including motion, RVT and HR did 
not appear to be major driving forces for the observed ensemble changes. The temporal 


changes did not correlate with global signal itself (Figure 2B). 


Evaluating Concordance among R-fMRI Indices: Voxel-wise Analyses 


From a static perspective, at each voxel, we calculated the Kendall's W (coefficient of 
concordance) of the 5 R-fMRI indices (fALFF/ReHo/VMHC/DC/GScorr) across the participants. 
As demonstrated in the concordance map (Figure 3A), the gray matter regions exhibit high 
concordance across individuals, including cerebral cortex and subcortical regions. In contrast, 
the white matter regions demonstrated low concordance. From a temporal dynamic perspective, 
we calculated the voxel-wise Kendall's W of the 5 R-fMRI indices across time windows within 
each individual. As demonstrated in the one-sample t-test on the Z-standardized within- 
individual concordance maps (Figure 3B), these analyses revealed that the voxel-wise 
concordance maps are highly similar to the static concordance. The rank order of the static 
concordance (Figure 3C) and the rank order of the dynamic concordance (Figure 3D) 


demonstrated a high correlation (r = 0.84, Figure 3E). 


Importantly, the voxel-wise concordance changes with age. As indicated by age effects on the 
concordance map, the brain showed widespread reductions in functional concordance as age 
increases, especially the subcortical areas (Figure 4A). To quantify the overall functional 
concordance, we averaged the dynamic voxel-wise concordance across gray matter voxels, and 
found this mean functional concordance index is negatively correlated with age across 


participants (r = -0.42, Figure 4B). 


Evaluating Spatial Concordance among R-fMRI Indices: Volume-wise Analysis 

Statically, we calculated the concordance (Kendall’s W) of the five measures across voxels for 
each individual, finding that volume-wise static concordance varies across participants, and is 
negatively correlated with age (r=-0.49, Figure 5A). Dynamically, the same concordance among 
the five R-fMRI indices was computed for each time window, yielding a concordance time-series 
for each participant (see Figure 5B for an example). As expected, the mean dynamic volume- 


wise concordance was highly correlated with the static concordance (r = 0.98, Figure 5C), and 


also negatively correlated with age (r=-0.49, Figure 5D). Of note, the CV of dynamic volume- 
wise concordance is positively correlated with age (r = 0.21, Figure 5E), suggesting that with 


age, the relationships among the various measures are slightly less stable. 


Understanding Low/High Concordance 

To understand the functional significance of concordance among different indices, we first 
divided participants into five groups based on volume-wise static concordance, and then 
compared the highest and lowest group using a two-sample t-test. Participants with the highest 
20% volume-wise static concordance demonstrated much higher values of all the five R-fMRI 
indices than those with the lowest 20% concordance (Figure 6A). This suggests functional 


concordance level also reflects traits differences in intrinsic brain activity properties. 


To further investigate relationships between concordance and the strength of R-fMRI indices, 
we carried out temporal dynamic analyses, which differentiated high and low concordance 
states on a within-individual basis. Specifically, for each participant, we sorted the windows 
based on their concordance score and we calculated the mean map for the highest 20% and 
lowest 20% of windows. For each R-fMRI index, this yielded two maps for each participant (high 
concordance, low concordance). Paired t-tests were employed to identify systematic changes in 
the R-fMRI indices across the two concordance states. The windows with high concordance 
demonstrated higher scores for all functional aspects (fALFF/ReHo/VMHC/DC/GSCorr) relative 
to those observed with low concordance (Figure 6B), suggesting high and low energy states 
existing in the brain. When looked at individually for each of the 7 large-scale networks 
previously defined by Yeo et al., the limbic network emerged as a notable exception, showing 


only modest increases in R-fMRI indices from low to high-energy states. 


Given the seemingly broad nature of increases in connectivity during high concordance states, 
we pursued one additional question — namely, whether increases in connectivity were specific to 
connections within network, or extend to those between networks? To answer this question, we 
investigated the intra-network functional connectivity and inter-network functional connectivity 
using the Yeo 7 networks (Yeo et al., 2011). For a pre-defined network, we calculated the mean 
time course for that network, and then correlated with all the voxels within another network. 
These correlations were averaged to represent the general correlation between the two 
networks. The correlation between the mean time course and all voxels within the same network 
can be averaged to represent the within-network functional connectivity. In general, within 
network functional connectivity is higher than between-network connectivity (Figure 7). For 
between-network connectivity, the limbic network demonstrated the lowest connectivity with 
others. The observation that both intra- and inter-network correlations were higher in the highest 
20% windows than the lowest 20% windows (Figure 7), suggests that both within- and between- 


network connectivity were increased during higher concordance states. 


DISCUSSION 

The present work provided a comprehensive examination of inter-individual variation and intra- 
individual temporal variation for commonly used measures of intrinsic brain function. We found 
that definitionally distinct R-fMRI indices of intrinsic brain function tend to exhibit a relatively high 
degree of co-variation within and between individuals. When taken as a measure of intrinsic 
brain function, inter-individual differences in concordance for R-fMRI indices appeared to be 
stable, and negatively related to age (i.e., functional concordance among indices decreases with 
age). In attempting to understand the functional significance of concordance, we noted that 
higher concordance was generally associated with higher strengths for the various R-fMRI 


indices. Finally, temporal dynamics analyses revealed that high concordance states are 


characterized by increased within- and between-network functional connectivity, suggesting 


more general variations in network segregation and integration. 


Our findings that differing R-fMRI indices exhibit commonalities in their patterns of variation 
across time and individuals have multiple implications. From a methodological perspective, they 
draw attention to questions regarding how to select an R-fMRI index for usage in a given study, 
as well as how to compare findings across studies that examine individual and group 
differences using different indices. The present work suggests that while differing indices share 
variation across individuals (and time), there was also unique variance associated with each of 
the indices — particularly in subcortical and limbic areas. As such, while it may remain justifiable 
for investigators to select one index or another based upon its congruence with a theoretical 
construct of interest, it is not fair to assume the differing measures are independent of one 
another. Perhaps more importantly, it is worth noting that we found that the concordance among 
the indices may in fact serve as a unique aspect of intrinsic brain function, which possesses 


stable trait differences among individuals. 


Beyond insights into study or analytic design, we believe the findings of the present work have 
potentially profound implications for our understanding of intrinsic brain function. In particular, 
the present work found that as concordance among measures varies over time, the degree of 
functional integration in the brain indexed by the various indices changes as well. This is 
analogous to the high efficiency states and low efficiency states defined by Zalesky et al. 
(Zalesky et al., 2014). For periods of high concordance, we found that notably different indices 
of connectivity (i.e., ReHo, DC, GCoRR, VMHC) and oscillatory power all showed maximal 
values, suggesting heightened functional integration. Consistent with this notion, connectivity 


both within and between large-scale functional networks was greater in these states. In contrast, 


low concordance periods were generally characterized by lower R-fMRI indices and decreased 


connectivity within and between networks. 


Examination of the correlation between the concordance among R-fMRI indices and 
physiological signals (e.g., cardiac, respiratory) failed to find any meaningful relationships, 
increasing suspicions of a potential neural etiology for this phenomena. Prior findings suggested 
that wide-spread fluctuations in fMRI signals (i.e., global signal) are correlated with spontaneous 
fluctuations detectable in the electrophysiological local field potentials (Scholvinck et al., 2010). 
Additionally, dynamic changes in functional connectivity were associated with variations in the 
amplitude of alpha and theta oscillations measured by electroencephalography (Chang et al., 
2013), and are linked to wakefulness (Tagliazucchi and Laufs, 2014). Consistent with these 
prior works, our findings suggested that global neural signals exist in the brain, and their 
spontaneous variations over time result in fluctuations in the connectedness of brain regions — 
both within and between functional systems, as well as locally (i.e., ReHo) and globally (i.e., DC, 
GSCorr). Additionally, our work suggests that these fluctuations in connectivity across brain 


regions are linked to fluctuations in the prominence of low frequency fluctuations. 


A key question is what drives inter-individual variation in concordance among R-fMRI indices. 
The present work draws attention to age as a contributing factor, finding progressive age-related 
decreases across the lifespan (age: 8-86). However, this only provides a clue into underlying 
etiology, the question that what age-related factor may be driving this relationship remains open. 
One possibility is the inter-individual variation in morphology (Yang et al., 2016), as gray matter 
volume changes with age (Sowell et al., 2003); such changes tend to show more complex 
patterns than the linear decreases with age observed in the present study. Second, is the 
possibility that it may reflect age-related changes in brain composition or neurochemistry. In this 


regard, it is important to note that iron, an essential precursor for dopamine synthesis, exhibits 


age-related increases in brain concentrations throughout the lifespan (Acosta-Cabronero et al., 
2016). A third consideration is recent work showing age-related decreases in brain temperature, 
which is thought to reflect by brain metabolism and cerebral blood flow (Sakai et al., 2011). 
Future work, integrating R-fMRI with other imaging modalities has the potential to differentiate 


between these possibilities, as well as to potentially reveal others. 


The present study revealed a temporal dynamic pattern of regional and global aspects of 
intrinsic brain activities. This dynamic nature is consistent with previous reports that focused on 
dynamic functional connectivity patterns (Chang and Glover, 2010; Sakoglu et al., 2010), based 
on which intrinsic brain activities could be decoded into several states (Allen et al., 2014; 
Leonardi et al., 2013); see reviews in (Calhoun et al., 2014; Hutchison et al., 2013; Tagliazucchi 
and Laufs, 2015). Interestingly, high-level cognitive regions, including the default mode network 
and dorsolateral prefrontal cortex, demonstrated that temporal variations were lower in regional 
activity measures (i.e., fALFF), but higher in functional connectivity indices (i.e., ReHo and DC). 
These regions are cortical hubs, and thus are expected to exhibit greater stability. This is 
supported by the high test-retest reliability of these regions (Zuo and Xing, 2014), and further 
supported by the lower temporal variability in regional strength. However, these regions are 
highly connected to other brain regions to process multimodal information (Buckner et al., 
2009); thus the connections should be flexible. This is supported by high temporal variability in 
functional synchrony measures in the current study, as well as by findings suggesting that these 
regions consistently forming dynamic functional connections (Zalesky et al., 2014). In contrast, 
the lower-level primary regions showed higher temporal variability in both regional strength and 
functional synchrony (high CV in fALFF, ReHo and DC), while with relatively lower actual 
regional activity (low mean in fALFF). This indicates the lower-level regions keep less active but 


more flexible states during resting-state. 


Several limitations in the present work merit consideration. First, the R-fMRI indices examined 
here were considered to reflect different functional aspects of intrinsic brain activity. There 
would be other aspects beyond the 5 R-fMRI indices examined in the present work, their 
variations may provide leverage for a deeper understanding of the physiological processes 
underlying the resting state signal. Second, we found the ensemble changes of intrinsic brain 
activity, which could be divided into high-integration and low-integration states. As simultaneous 
electrophysiological recordings were not available for these datasets, we can only speculate 
how the neuronal ensemble firing related to these fluctuations in functional concordance. Third, 
we found how functional concordance is related to age in healthy participants, while its 
generalizability to clinical populations remains unknown. Further studies should apply these 


methods to neuropsychiatric disorders to test their sensitivity and specificity in brain disorders. 


In sum, the present work provided a comprehensive examination of the functional concordance 
of intrinsic brain indices from inter-individual variation and temporal dynamics perspectives. The 
current study draw attention to questions regarding how to select an R-fMRI index for usage in a 
given study, as well as how to compare findings across studies that examine differences 
between individuals or groups using different indices. Our work suggests global neural signals 
exist in the brain, and their spontaneous variations over time result in fluctuations in the 


connectedness of brain regions. 
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FIGURE LEGENDS 

Figure 1. The static and temporal dynamic patterns of R-fMRI indices. T scores of static (A), 
temporal mean of dynamic (B), and temporal coefficients of variation of dynamic (CV, C) R-fMRI 
indices. Average t values of Z-standardized temporal CV of dynamic R-fMRI indices were 
computed within each of the seven networks defined by Yeo et al (2011) and plotted using radar 
plots (D). SomMot: Somatomotor; DorsAttn: Dorsal Attention; VentAttn: Ventral Attention; 


Control: Frontoparietal Control. 


Figure 2. Global-level concordance among R-fMRI indices. (A) Static concordance: Pearson's 
correlations between pairs of global mean of static R-fMRI indices were computed across 
participants. (B) Dynamic concordance: Pearson’s correlations were computed between pairs of 
global mean time series derived from dynamic R-fMRI indices within a participant and then 
averaged across participants. For dynamic concordance, the contributions of nuisance variables 
were also examined: motion, global signal (GS), respiration volume per time (RVT) and heart 


rate (HR). 


Figure 3. Functional concordance among five R-fMRI indices: voxel-wise analyses. (A) Static 
concordance among R-fMRI indices computed across participants (Kendall’s W); (B) Dynamic 
concordance among R-fMRI indices (T scores) computed in two steps: first, in each participant, 
Kendall’s W was computed across time window among five R-fMRI indices; then one-sample T- 
tests were performed to collapse the concordance maps (after Z-standardization) across 


participants; (C) Rank order of A; (D) Rank order of B; (E) Spatial correlation between C and D. 


Figure 4. Age effects on dynamic voxel-wise within-individual concordance (Kendall's W). (A) 
Wide-spread reductions in functional concordance as age increases. Of note, the result is 


corrected at voxel-level P < 0.001 (T > 3.1 given the degree of freedom = 169) and cluster-level 


P < 0.05. (B) Mean functional concordance index (dynamic voxel-wise concordance averaged 


across gray matter voxels) is negatively correlated with age across subjects. 


Figure 5. Functional concordance among five R-fMRI Indices: volume-wise analyses. (A) Static 
volume-wise spatial concordance among R-fMRI indices was plotted as a functional of age; (B) 
An exemplar dynamic volume-wise concordance among R-fMRI indices for a given subject (one 
concordance time-series/participant); (C) Correlations between static and mean dynamic 

(averaged time windows) volume-wise spatial concordance; the mean (D) and the coefficient of 


variation (CV: E) of dynamic volume-wise spatial concordance was plotted as a function of age. 


Figure 6. Differences in magnitude of R-fMRI indices between high and low concordance 
participants or time windows (A) T cores obtained via comparing the static R-fMRI indices 
between participants with highest and lowest 20% of concordance (static volume-wise spatial 
concordance) were plotted for each voxel on surface map (top row). Of note, the result is 
corrected at voxel-level P < 0.001 (T > 3.2 given the degree of freedom = 68) and cluster-level P 
< 0.05. The average t scores of each of the seven Yeo et al., (2011) networks were plotted in 
radar chart (bottom row). (B) For each participant, the mean R-fMRI index maps for the windows 
with highest 20% and lowest 20% concordance (dynamic volume-wise spatial concordance) 
were calculated, yielding two maps (high concordance, low concordance), and then paired t- 
tests were performed between high and low concordance maps across participants. The result 
is corrected at voxel-level P < 0.001 (T > 3.1 given the degree of freedom = 172) and cluster- 
level P < 0.05. SomMot: Somatomotor; DorsAttn: Dorsal Attention; VentAttn: Ventral Attention; 


Control: Frontoparietal Control. 


Figure 7. Intra- and inter-network functional connectivity between high and low concordance 


state. (A) Mean functional connectivity averaged among the highest 20% windows of dynamic 


volume-wise spatial concordance. (B) Mean functional connectivity averaged among the lowest 
20% windows of dynamic volume-wise spatial concordance. (C) Paired T-tests on functional 


connectivity between the highest and the lowest 20% windows. 
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